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ABSTRACT 



The complete two-loop correction to the quark propagator, consisting of the 
spider, rainbow, gluon bubble and quark bubble diagrams, is evaluated in the non- 
covariant light-cone gauge (leg), n-A a (x) = 0, n 2 = 0. (The overlapping self-energy 
diagram had already been computed.) The chief technical tools include the powerful 
matrix integration technique, the n*-prescription for the spurious poles of (q-n)^ 1 , 
and the detailed analysis of the boundary singularities in five- and six-dimensional 
parameter space. It is shown that the total divergent contribution to the two-loop 
correction S 2 contains both covariant and noncovariant components, and is a lo- 
cal function of the external momentum p, even off the mass-shell, as all nonlocal 
divergent terms cancel exactly. Consequently, both the quark mass and field renor- 
malizations are local. The structure of £2 implies a quark mass counterterm of the 
form Sm(lcg) = ma s C F (S + a s W) + 0(5f), a s = g 2 r(e)(47r) e ~ 2 , with W depending 
only on the dimensional regulator e, and on the numbers of colors and flavors. It 
turns out that 5m(lcg) is identical to the mass counterterm in the general linear 
covariant gauge. Our results are in agreement with the Bassetto-Dalbosco-Soldati 
renormalization scheme. 
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1 Introduction 



A quarter of a century ago the light-cone gauge was a gauge "to fortune and to fame 
unknown" . It was regarded as a freakish member of the family of axial-type gauges 
that existed more by accident than by inventive planning [|J. Today, the light- 
cone gauge enjoys respect and a privileged status among noncovariant gauges for 
the following two reasons: first, computations in the light-cone gauge have proved 
to be meaningful, both at one and two loops. Second, the renormalizability of 
light-cone QCD, as demonstrated in this article, has finally been established at 
two-loop order by explicit computation. Specifically, we shall discuss here in some 
detail the renormalization of the quark propagator to two loops in the light-cone 
gauge n- A a = 0, n 2 = 0, where A® denotes the gauge field and = (no, ft) is an 
arbitrary, but fixed, four- vector, // = 0, 1, 2, 3 0. The diagrams for this process are 
depicted in Figure 1. We note that the results for the one- loop quark self-energy 
function (Fig. If) and for the overlapping self-energy function (Fig. lb) have already 
been reported in the literature 0, ||. These two diagrams are included here for 
completeness. 

Whereas computation of Fig. If was straightforward, at least in hindsight, eval- 
uation of the overlapping quark self-energy |4], |5| required the introduction of a new 
procedure, called the matrix integration technique. We recall that in this procedure, 
the two momentum integrals of J d 2uJ q J d 2u} kf(q, k), where 2uj denotes the complex 
dimensionality of space-time, are integrated over 4o;-dimensional space in a sin- 
gle operation. We further recall that the biggest advantage of the matrix method 
is the ability to execute the momentum integrations exactly and in closed form. 
With the momentum integrals conveniently "out of the way" , we can then concen- 
trate on the wide variety of new parameter singularities which is so characteristic 
of noncovariant-gauge multi-loop integrals. It turns out that the matrix method 
enables us to handle these new parameter singularities in a consistent and unam- 
biguous manner. By comparison, multi-loop integrals with fewer and less severe 
parameter singularities may, in general, be evaluated by means of the nested method 
H> H> Hi 01 • I n this traditional approach, the 2k>momentum integrations are carried 
out sequentially. 

We finally observe that the matrix integration technique works for covariant 
and noncovariant gauges alike, and regardless whether the integrals are massive or 
massless. We shall have occasion later in this article to highlight specific technical 
features of this powerful technique. 
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e) quark-bubble f) one-loop 

Figure 1: One- and two-loop quark self-energy diagrams for QCD. Wavy lines denote 
gluons, straight lines denote quarks, p, q, and k are four-momenta. 
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It is well known || || that in covariant gauges, QCD may be renormalized by a 
redefinition of the masses, coupling constants, and field normalizations. However, for 
noncovariant gauges, such as the light-cone gauge, the situation is more complicated: 
in addition to the same types of counterterms that arise in covariant gauges, for 
light-cone QCD we also require n^-dependent counterterms. These noncovariant 
counterterms are different from any of the terms in the original Lagrangian density. 
On the other hand, and in contrast to the covariant case, no physical measurements 
are needed to fix the finite parts of these new noncovariant counterterms. Instead, 
the finite parts are constrained by the requirement that all physical observables 
be Lorentz-covariant. In a sense, the presumption of covariance of the observables 
constitutes an infinite set of physical constraints on the noncovariant counterterms. 

In 1987, Bassetto, Dalbosco, and Soldati (BDS) jR], [□]] proposed a renormal- 
ization scheme for light-cone QCD in which the particular types of noncovariant 
counterterms - permitted by gauge symmetry and Lorentz symmetry - are absorbed 
into the original Lagrangian density by means of noncovariant renormalizations of 
the quark and gluon fields. Unlike the covariant-gauge case, where the renormaliza- 
tion factors are scalars, in the BDS scheme the renormalization factors are matrices. 
The result is that the different spin- components of the quark and gluon fields are 
renormalized differently. The present two-loop calculations are a clear vindication 
of the BDS renormalization scheme. 

The plan of paper (II) is as follows. In Section 2 we derive the integrands of the 
light-cone integrals for the various diagrams of Figure 1, and then break them down 
into approximately 80 simpler integrands. In Section 3 we illustrate, by means of 
examples, methods for handling the various integrations. The explicit values for 
the divergent parts of all integrals are tabulated in the Appendix. These results 
are analytic and allow for general quark masses. The counterterms needed for the 
renormalization of the quark self-energy to two loops in the light-cone gauge are 
derived in Section 4. It turns out that all required counterterms are local (that 
is, polynomial in the external momentum), and that their coefficients satisfy the 
relationships implied by the BDS renormalization scheme. The paper concludes in 
Section 5 with a summary of our results and their significance. 



2 Derivation of the Integrals 

The Lagrangian density for light-cone QCD reads 

£ = -\F;„F a ^ + ^ a (i$ - m)i, a + gT^ a 4 a ^ + Ax + £ gh + At , (1) 
where external source terms and quark flavor indices have been suppressed, and 
F% = d^A* - d v A% + gf ahc A\A% , = . . 3, d, = d/dx», 
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A a Ax) — gluon fields, a,b,c= gluon indices = 1 .. 8 for SU(3), 

"0a( x ) = quark fields, = color indices = 1,2,3, ip = V^7°j 

jfi = j- A, $ = •y-d, etc., 7 M = Dirac matrices, (7 ) 2 = 1, 

m = quark rest mass, g = coupling constant, 

T£g = generators of SU(3), normalized such that [T a ,T b ] = if abc T c , 

jabc _ an tisymmetric structure constants of SU(3) (with f abc f bcd = 35 ad ), 

^fix = gauge-fixing term = — (n- A) 2 /(2X), with A to be taken to zero later, 

£ g h = ghost terms = —r] a (5 ab n-d — gf abc n-A c )r] b , if 1 = ghost fields, and 

£ ct = counterterms. 

Repeated indices imply summation, and we use a +,—,—,— metric for Minkowski 
space. 

2.1 The Feynman Rules 

The Lagrangian (p]) leads to the following Feynman rules for light-cone QCD |TT 
0,|: 



Quark propagator: 

tfaf,S(p), With S(p) = J + ™-n > ^>0, (2) 

where p M is the quark 4-momentum, a and (3 are color indices as in (|TJ) , and the i6 
term is Feynman's prescription for avoiding a singularity when p 2 = m 2 . (We let 9 
go to zero after Wick rotation.) 

Gluon propagator: i5 ab G ^(q), with 

BM - -iU L - + , e > 0, (3) 

q z + tV \ n-q [n-q) z J 

where g M is the gluon 4-momentum. The gauge-fixing parameter A is now taken to 
zero, causing the third term in parentheses to drop out. 

Quark-quark-gluon vertex factor: igj^T^. 

3-gluon vertex factor: gf abc [(p — q) p 9^ u + (<? ~ k)^g up + (k — p) I/ <7 w> ], where 
p, q, k are the incoming 4-momenta of the attached gluons. 
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4-gluon vertex factor: -ig 2 [f abe f cde ( g w g ^ -gf» g »p)+f<«* fde^nv gP a _ gi w g u P ^ + 

jade jcbe^g^pgua _ g/j,v g(xj^ 

Ghost propagator: 5 ab /(n-k), where is the ghost 4-momentum. 

Ghost-ghost-gluon vertex factor: —igf abc n fl , where \x and a match indices of 
the attached gluon, while b and c match indices of the attached ghosts. Note that 
ghosts "decouple" in the light-cone gauge, because the ghost-ghost-gluon vertex 
factor is orthogonal to the gluon propagator (H) when A = [14], |J . 

Counterterm vertex factors: to be determined in Section 4. 



To construct the dimensionally regularized Green function for a diagram, we 
first impose conservation of momentum at each vertex, so that, apart from external 
momenta, there is only one independent momentum per loop. (In Figure 1, these 
momenta are denoted by q and k.) We then form the product of the propagators 
and vertex factors for all internal lines and vertices, divide by (2ir) 2uJ for each loop, 
and integrate the resulting expression over the 2o;-dimensional space of each loop 
momentum. Finally, we multiply the integral by —1 for each internal quark loop, 
divide by r! for each pair of vertices connected by r gluon lines, and also divide by the 
number of permutations of the vertices which leave the diagram invariant (for fixed 
external lines) [11, |I5 |. In Figure 1, only diagram e (quark-bubble) has an internal 



quark loop, only diagram d (gluon-bubble) has a pair of vertices connected by more 
than one gluon line, and no diagram is invariant under a non-trivial permutation of 
its vertices. 

The procedure just described gives only the "amputated" Green function, since 
it includes no factors for the external lines of the diagram. 



2.2 The Integrals and the n*-Prescription 

To illustrate the application of the above Feynman rules, we may immediately write 
down the amputated Green function for the one-loop diagram of Figure If: 

SaJliip) = ~f 2 ^2? K j M ^ S (P ~ G M (4) 

where J M denotes integration over Minkowski space, and S and G are the propaga- 
tors defined in equations (0) and @ (with A = 0, and m equal to the rest mass of 
the external quark). To avoid a singularity in G when n-q = (in this and other 
integrals), we shall use the n*-prescription (ML-prescription) [16], in which 

-1- -> lim W *- g - =hm 2n *-\. a , 9>0, (5) 
n-q o^o n*-qn-q + iv e^o n*- n q^ + 2tt> 
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where n* is a new light-like 4- vector (n* 2 = 0), with n*-n > in Minkowski space, 
and 

n*-qn p + n-qn* 

9||a» = > <?±m = 9/x ~ I 6 ) 

for any 4-vector g M . The same prescription applies to (n-k)' 1 , with the same n*, 
but with k in place of q. Prescription @ was subsequently recovered in the context 
of canonical quantization by Bassetto et al. [0. Unlike the old "principal value" 
prescription, the n*-prescription is consistent with both Wick rotation and power- 
counting 0,0, 0,0. 



The amputated Green functions for the five two-loop diagrams of Figure 1 possess 
the following structure: 

Rainbow diagram: ig 4 T^T^ s T s a K (27T)-^I a , with 

I a = / d 2 "q [ d^k^S(p-q)YS{p-q-k)YS{p-qh a G^{q)G up (k). 
Jm Jm 

Overlapping diagram: ig 4 T®pTp J T" s Tg K (2 / 7r)~ 4u ' 7 b , with 

J b = / d 2 "q [ d 2 "kYS(p-q)YS(p-q-k)YS(p-k)YG, p (q)G ua (k). 
Jm Jm 

Spider diagram: g* f abc T^T b Pl T^{2-K)-^ I c , with 

d 2 "q [ d^k^S(p - g)7^(p - k)Y G^q) G vv {k - q) G pT (k) 
m Jm 

■[(k - 2g)V" + (q ~ 1kfg vp + (q + k) v g^\. 

Gluon-bubble: \ig 4 f abc f bcd T* p T^{2n)-^ I Al with 

h=f d 2u >q [ d 2 "k^S(p-q)^G^(q)G uv (k-q)G pT (k)G^(q) 
Jm Jm 

.p _ 2q)Pg^ + (q- 2k)»j"> + (q + k) VI 
•[(& - 2g) T ^ + (q - 2kYg r,T + (q + 

Quark-bubble: -tg 4 T^ 5 T^T b K (27i)-^I e , with 7 e = 

(7) 

where is the rest mass of the quark in the inner loop of Figure le. 
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The first step in the evaluation of integrals I a to Id is to substitute for S and G 
from equations (|2|) and (|3]) , and expand the numerator of each integrand into a sum 
of products. In order to avoid terms whose integrals diverge as q — > 0, the quark- 
bubble integral I e will be handled differently from the other integrals, as discussed 
in Section 3.2. 

Before integrating, we simplify the integrals by applying Wick rotations. This 
procedure is valid because both Feynman's id prescription as well as the n*-prescription 
lead to poles in the second and fourth quadrants only. In this way, the integral 
for instance, becomes 



(2tt) 2 ^ Je [(p- q) 2 + m 2 }q 2 \ ^ n-q J' 

where //, v = 1 ... 4, q = iq±, p = ip±, n = in±, 7 = ry 4 , and J E denotes 
integration over the Euclidean space spanned by q± ... q±. 



2.3 Expansion, Reduction, Transformation, and Tadpoles 

From the Feynman rules, we see that Green functions in general are integrals of 
rational functions of the loop momenta. For integrals 7 a to 1^ above, it would 
appear that some terms in the numerators of the integrands could have degrees as 
high as 11 in q and k together (after up to four applications of the n* -prescription 
d^D). Fortunately, however, we may reduce the maximum degree to five, and in most 
cases to three, by means of cancellations between terms within each integral, and 
between numerator and denominator factors within many of the terms. Because of 
the large number of terms involved, some computer assistance is advantageous. 

Expanding the numerator of each integrand into a sum of products, and then 
following the procedure outlined in paper (I), we obtain a list of well over 100 distinct 
terms to be integrated. Fortunately, we can shorten this list still further by applying 
transformations, such as k — > q — k, q — > q + k, and/or q <-> k, to selected terms. 
The transformation q — > q + k, for instance, yields 

(n-q — n-k)(q — k) 2 [(p — q) 2 + m 2 ]k 2 n-qq 2 [(p — q — k) 2 + m 2 ]k 2 ' 

where the expression on the right just happens to be already on the list of terms to 
be integrated. Notice that the factor (n-q — n-k) can thus be eliminated from all 
denominators. 

Having carried out these various cancellations and transformations, we find that 
some of the integrands factor into separate q and k dependent parts. In some cases, 
one of these two one-loop integrals vanishes, since its integrand is antisymmetric 
under q — > — q or k — >• —k. In many other cases, one of the one-loop integrals 
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corresponds to a massless tadpole, and likewise vanishes in the context of dimensional 
regularization [ITSL O . 



2.4 Power Counting 

For renormalization, we require only the divergent parts of integrals J a to J e , so we 
shall drop all terms whose integrals can be shown by power counting to converge 
when uj = 2. For a covariant gauge, Weinberg's theorem |2(| [15|] tells us that a Feyn- 



man integral is UV- convergent if its integrand, including the measure d 4 qd A k . . ., is 
of negative degree with respect to every non-empty subset of the loop momenta 

q,k, We have the same rule for the light-cone gauge, except that we must also 

consider subsets that include only the "transverse" part - the part orthogonal to n 



- of one or more of the loop momenta [[Ti l. 

To illustrate power counting, consider the Euclidean-space integral 

ifin-pd 2uJ qd 2uJ k 



n-q[(p — q) 2 + m 2 ](q — k) 2 (n-q — n-k)[(p — k) 2 + m 2 ] ' 

arising from the spider diagram (Figure lc). When uj = 2, the integrand has degree 
—2 in q, —1 in k, and in q and k combined. Hence, by the rule given above, 
this integral may be divergent. As it turns out, however, this particular integral is 
actually convergent, because the leading-order part of the integrand for large \q\ and 
| A; | is antisymmetric under q% <-> k% «-> k± (we have taken ri\ = n 2 = n\ = n* 2 = 
for simplicity, so that q^ = (0,0,^3,^4) and k\\^ = (0, 0, k$, k±)). There are other 
such "borderline" integrals which are convergent for the same reason. 

In the Appendix we have summarized the divergent terms which remain to be 
integrated, for each of the integrals I a to I e . We have also listed there the various 
integrated divergent parts for each individual term. In the next section we shall 
demonstrate how these results were obtained. 



3 Integration Methods 

Several different approaches to the evaluation of multi-loop Feynman integrals ap- 
pear in the literature P, [7], |2T], E3|; they will not be reviewed here. In most 
methods, one begins by writing the denominator of the integrand as an integral, 
using a parametrization formula. We shall use the formula known as Schwinger's 
representation: 

1 f°° c^dcti f°° al^da 



Jo -TfuWo "TfVT " ' 6XP ^ ~ - « 2 F 2 . . . ). (9) 
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The factors Fi,F 2 , . . . must have positive Real parts; hence, before using formula 
(|D, we apply the n*- prescription to any factors of n-q and n-k in the denominator, 
perform Wick rotations, set 9 to zero, and take the factor(s) n*-n from Eq. (||]) outside 
of the integral. The component P4 of the external momentum is regarded as Real 
until after integration has been completed. 



3.1 The Nested Method 

After parametrization, it is necessary to integrate over both the loop momenta and 
the parameters. The question of the most suitable order for these integrations 
naturally arises. To take advantage of known one-loop results, we might try the 
nested method, in which we integrate first over one of the loop momenta (k, say), 
then over the parameters associated with factors involving k, then over the other 
momentum q, and finally over the remaining parameters. 

As an example, consider the divergent Euclidean-space integral 

h= f f jd^qd^k 

J J n-q[(p — q) 2 + m 2 ][(p — q — k) 2 + m?]k 2 ' 

arising from the overlapping and rainbow diagrams. Applying formula @ to the 
A;-dependent factors only, we obtain 

r tk d 2u ' q r 00 f°° r / 

h = / -± / dai / da 2 / d 2w k exp ( - a-Mp - q - k) 2 + m 2 } - a 2 k' 

J n-qQ Jo Jo J v 

where Q = (p — q) 2 + m 2 . Next we change variables from ai,a 2 to A = ai + a 2 . 
x = a.i/ A, and complete the square in the exponent to get 

h = [t^-^ ( l dx [°°AdA [d 2uJ kexp(-AUk-xp + xq) 2 + H]), (11) 
J n-qQ Jo Jo J v ' 



with H = x 



(1 — x){p — q) + m . The k and A integrations may then be carried 



out with the help of the well-known Gaussian and Gamma integrals: 

J exp {-Ar 2 ) d 2uJ r = Cjj" , A > 0, (12) 

and 

r°° V(d 4- 1 1 

A d e -AH dA = g_r_ i ; H>0 ^ Rerf> _i ; (13) 



respectively, with r Ph = k ll — x(p — g) M in this case. Accordingly, we obtain 

I x = 7r"r(e) fi^l f'dxH-t; e = 2-uj. (14) 

J n-aUJo 
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From power counting, we expect the k integral in Eq. fllPl) to be well denned only 
if u < 2. In fact, we required just this condition in order to complete the integration 
fll3|) which produced the divergent factor T(e) in Eq. (0). The q integral in Eq. ( |i~i|) 
also diverges as e — > 0, according to power counting, so altogether we expect I\ 
to have a double pole at e = 0. This expectation will be confirmed by explicit 
calculation. 

Before we can complete the integration, we must decide how to deal with the 
g-dependent factor H~ e in Eq. (|H]) . Since we expect I\ to have a double pole at 
e = 0, and since we are only interested in finding the divergent parts, we might try 
integrating only terms up to order e from the exponential series 



H 



1-elnif + Me\nH) 2 - ... 



1 — e In x — e In 



[1 - x)(p - q) 2 + m 2 } + |(elnx) 2 + . . . 



x 



-e 



(1 — x){p — q) 2 + mf 



l + 0(e 2 ). 



(15) 
(16) 

(17) 



If we could drop the 0(e 2 ) terms, we could immediately integrate over x, and then 
complete the q integration from Eq. flT3|). Unfortunately, the series (|I5) cannot be 
integrated term-by-term, since it does not converge uniformly with respect to q, as 
explained in paper (I). 

The convergence of series ([T^) is non-uniform partly because H goes to infinity 
as \q\ goes to infinity. One way of solving this problem is to factor out the large 
\q\ behaviour before using a series, provided we do so without either creating new 
convergence problems, or generating terms that we cannot integrate. In the current 
example, we can extract a factor of Q from H before using the exponential series, 
so that Eqs. ([151 ) to ([H]) become 



H 



Q- 



eln 



H 

Q 



eln 



H 

Q 



Q 1 — elnx — eln 



^1 — x){p — q) + m 

Q 



+ -(elnxr + • • • 



^1 — x)(p — q) + m 

' Q 



1 + 0(e 2 ) 



(19) 



with Q = (p- 

exceptnea.1 x 



■q) 2 + m 2 . Convergence of this new series remains uniform as \q\ — > oo, 
= and x = 1. Fortunately, however, these remaining non- uniformities 
cause no trouble, provided we integrate Eq. ([19|) only. At x — 0, the x~ e term 
in this equation gives the correct contribution to the pole parts of I±, while the 
contributions from the second and third terms in parentheses cancel. Similarly, as 
x — > 1, |g| — > oo, the second term in parentheses gives the correct contribution, 
while the contributions from the first and third terms cancel. (One may check these 
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claims by using expansions in powers of x near x = 0, and 1 — x near x = 1.) Hence, 
the 0(e 2 ) term in Eq. (|19|) does not contribute to the pole parts of I\. 



Integrating Eq. ([19|) over x, and substituting into Eq. (|T4|), we find that 



h 



7t u T(e) f^d 2uJ q 
n-qQ 



1 + e m 2 [Q- e - (m 2 )' 



Q ( 



(p - <i¥ 



+ finite, 



(20) 



where "finite" refers to terms which do not diverge as e — > 0. By power counting, 
we see that the q integration in Eq. ( p0|) diverges only for the first term in square 
brackets. Because of the divergent factor T(e) in front, we need both the divergent 
and finite parts of the q integral, but not parts of order e. Hence, we may set e = 
in the convergent second term, causing this term to vanish (fortuitously). Applying 
the n* -prescription (|5]), along with the parametrization formula (|9|), we obtain 



2?i7r a T(e)(l + e) f°° , f°° a%da 2 

h = ; ; — / dcti / — 

n*-n(l-e) Jo Jo r(l + e) 



n -qe 



- T d 2u) q + finite, 



2p"{l + e) 



^4 [\ e dy rA 1+e dA [n*-(r + y Pll +p ± )e- T d 2ul r + finite, (21) 
1 — e) Jo Jo J 



where A = a x + a 2 , y = a 2 /A, r^^q^- ypy - p ± ^, and 



T = aiiq 2 + a 2 Q 



A 
A 



Q\\ + vil ~ 2 yp-i + y(p 2 + m2 
r \ + y r ± + 2/(1 — y)p\ + y m ' 



The 4- vectors py, p±^, and so on, are defined in Eqs. @. Note that n*-\ 
P'Q — P\\'Q\\ +P±"9±) ip\\)± — 0, etc., for any 4-vectors p^ and q^. 



n -p, 



Since T is an even function of r, we may drop the term with n*-r from the 
integrand of Eq. (plf) . We also drop n*-p±, which is equal to zero, and factor the r 
integral into independent r\\ and r± parts: 



yn -pii e d r 



y n*-p exp ( — Ay (1 — y)p\ + ?ti 
x y exp (— Ar 2 ) d 2 ry / exp {-Ayr 2 ± ) d 2w ' 2 r ± . 

The integrations over r, A, and ?/ are then easily completed with the help of formulas 
(1121) and (ITT 



3.2 The Quark-Bubble Integral 7 e 



Let us try the nested method on the integral in Eq. (0). Since the inner loop of 
Figure le involves no gluons, the integral over k of the fc-dependent factors in (0) is 
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just the standard covariant-gauge result |15| : 
,awi Trace[7 CT (^ + m L )7^(^-^ + m L ) 



d^k 



m [k 2 — m\ + i9] \{k — q) 2 — m| + i6\ 

8i7r"T(e)(gV P - </V) f (x - x 2 )\m\ - (x - x 2 )q 2 ]~ e dx, (22) 

Jo 

with 9 > and e = 2 — u; as before. The above equation is exact: no terms of 
order e , or higher, have been omitted. (To verify this claim, one may use the trace 
theorems Jl5| to obtain 

Trace [y a {ty + m L )Y{ty - <t + m L )} = Ag ap (m 2 L - k 2 - q-k) + 8k a k p - Aq a k p - Ak a q p , 
and then carry out the k integration using some of the formulas and methods from 
Subsection 3.1.) 



Substituting the right-hand sides of Eqs. (|22| ) and (||) into Eq. (0), we get 

Ie = 8m - r(e ) f d 2^ 7 M (# zA± m ^ D f q ) f 1 ( x ~ dx (23) 
Jm (p — q) 2 — m? + %Q ^ Jo \m\ — (x — x 2 )q 2 ] e ' 

where D^i^q) = G pa (q)(q 2 g crp — q° q p )G pu (q) , and m is the rest mass of the external 
quark in Figure le. It follows from Eq. (|) that D^q) — > —G^q) as 9, A — > 
0, so that I e is proportional to the one- loop integral (||), except for the extra x 
integral within the q integral. We note that there is only one factor of q 2 now in the 
denominator (in D piJ (q)), thereby preventing the integral from diverging as q — > 0. 



After Wick rotation, Eq. (p^) becomes 

'.= 87r«T(e) [d 2 "q Ti ! " . J <\„. ) I " ' «!■/•■ 



Yty-j- ™)l u f 5 n^ + q^nA f 1 (x - x 2 ) 1 ( 
[(p — q) 2 + m 2 ] q 2 \ ^ n-q J Jo L e 

(24) 

with L = q 2 + m 2 L (x — x 2 ) . As in the example of Eq. (|10|), power counting tells 
us that integral (^) is well defined only if e > 0. Under this condition, {q 2 )~ e is a 
decreasing, concave- up function of q 2 (for q 2 > 0), from which it follows that 



d 



em, 



< (qT'-L-* < (^-L)^[(^)" e J = ^T^ ) ( 25 ) 

Using this inequality, we can show by power counting that the difference between 
integral (|24]), and the same integral with L replaced by q 2 , remains finite as e -> 
0. Since we are merely interested in the divergent parts of I e , we may make this 
replacement and pull the factor {q 2 )~ e out of the x integral. Integrating over x with 
the help of the formula 

^ T(c+l)T(d+ 1) 

/o ~ v ~ ~' ~~ T(c + d + 2 



x c {l~x) d dx = \T t \ . Rec,Red>-l, (26) 



and completing the g integration, we get the desired result. Notice that the divergent 
parts of I e will be independent of the mass of the quark in the inner loop, as m^ 
enters Eq. (|24]) only by way of L. 
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3.3 Review of the Matrix Method 



As seen in the preceding examples, in the nested method one tries to modify the inte- 
grand between the first and second momentum integrations, in a way that simplifies 
the final integration without changing the divergent parts of the result. Particular 
care must be taken with regard to the behaviour of the integrand near boundaries 
where the final integration diverges, such as \q\ —>■ oo, x — > in the examples. The 
greater the degree of divergence, the more closely the "simplified" integrand must 
match the exact one. This requirement becomes more and more challenging as the 
number of a parameters increases. 



For integrals such as 

fa tf^q^k 



e je n-q[{p — q) 2 + m 2 ] [(p — q — k) 2 + m 2 ] [(p — k) 2 + m 2 ] k 2 



(27) 



in which both q and k appear in more than two denominator factors each, it is 
easier to complete all momentum integrations before doing any parameter integra- 
tions. This approach is called the matrix method [|], |5|. Momentum integrations 
are straightforward with this method, because the combined 4u;-dimensional mo- 
mentum integral can be expressed as a derivative of a product of one-dimensional 
Gaussian integrals. In view of the importance of the matrix method for multi-loop 
integrals, we shall briefly review its main features. 

After Wick rotation and application of the n*-prescription, a two-loop light- 
cone integral takes the form of an integral over q and k of a polynomial P(q,k), 

say, divided by some non- negative quadratic factors F±, F 2 , Application of the 

parametrization formula @ then yields 



d 2» q l^ k EML = I da l d a 2 ... J[P(q,k)}; (28) 

JE r\t 2 JO JO 

J[P(q, k)} = Jd^q J d 2w k P(q, k) exp ( - ai F x - a 2 F 2 ...). (29) 



00 



Since the exponent is quadratic in q and k, we can rewrite Eq. (|29|) in the form 

J[P(q,k)] = J E d 4ul zP{q,k) exp(-zMz T + 2B-z-C), (30) 

where M is a Real 4o;x4o; matrix, B and z are Real 4co>-vectors, z = (^4,^4, k 3 ,q 3 , 
. . .), T denotes transpose, and the components of M, B, and C are functions of 
«i, «2 • • -j but not of q or k. 

For any given set of F factors, one may obtain explicit expressions for the com- 
ponents of M, B, and C by expressing the exponents from Eqs. ( p9| ) and ( |30| ) in terms 
of the components of q and k, and then equating corresponding coefficients. For the 
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two-loop integrals of the Appendix, with q\\ = (0, 0, q 3 , q^) and k\\ = (0, 0, k 3 , k±) as 
before, we find that 



M = A 



a tG 
tG 1-a 



a tG 
tG 1-a 



h tG 
tG A 



B 



A 



bp A 
bp 2 

PP2 



(31) 



where the dots in M denote 2u — 3 repetitions of the third 2x2 sub-matrix, and 
A(l — a), AX, A/3, AG, Ab, Ah, and Aa are the sums of the a parameters whose 
corresponding F factors include gij, qj_, —2p-q, 2tq-k, —2p-k, kj_, kp 
respectively. For example, if we label the denominator factors in integral (p7|) as F\ 
to F 5 from left to right (after application of prescription (|5|), with 2n*-q/n*-n taken 
into P), we have for this integral 



1 

h 



a - 
■ a 



(«i + a 2 + a 3 )/A, 
: (a 3 + «4 + a 5 )/A, 



\ = [3 = (a 2 + a 3 )/A, t = 1, 
b = (a 3 + a 4 )/A, G = a 3 /A, 



C = (a 2 + a 3 + a 4 )(p 2 + m 2 ) = A(b + X - G)(p 2 + m 2 ). 
The equations for a and 1 — a tell us that A = a± + a 2 + 2a 3 + 04 + 0:5. 



(32) 
(33) 



Since the matrix M is to be multiplied by z on both sides, it may always be 
constructed symmetrically. Accordingly, there exists a matrix R such that RMR T is 
diagonal and R T R = 1. Defining a new vector y so that z = yR + BM^ 1 , and taking 
P = 1, we find that Eq. (|30[) becomes 

J[l] = exp(BM- 1 B T -C)^|detR|d 4w y exp ( - yRMR T y T ) • (34) 

Since the integral in this equation is just a product of 4a; one-dimensional Gaussian 
integrals, we may use formula ([l"2j), along with |det R| = 1 and ni(RMR T )jj = 
det(RMR T ) = det M, to get 

J[l] = 7r 2 "(detM)" 1/2 exp(BM- 1 B T -C). (35) 

Notice that we never actually need to construct R. To find J[P] for general P(q, k), 
we differentiate Eq. ( pO] ) partially with respect to Bj to obtain dJ[P]/dBi = 2J[zjP], 
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and then apply this formula repeatedly to Eq. (|3~5"D to get 



J[zi] = TiJ[l] 

J[ZiZjZ k ] = 



TiT,- + (iM 



y 

H/l-i 



l^ + T^iM-^ + TfcdM- 1 ) 



(36) 



where T = BM 1 = (r 4 , s 4 , r 3 , S3, ... 



(37) 



Since the momentum integral J[P] is a linear functional of P, and z = (£4, g 4 , £3, g 3 , 
. . .), we can use the above equations to find J[P] for any polynomial P(q, k). 



In order to derive the two-loop integrals summarized in the Appendix, we sub- 
stitute from Eqs. ([31]) into Eqs. ( |35"D to fl37|) to obtain the following relations: 



J[%] = 8„J[1] 

J[n*-qn-k k^] = 



J[n*-qk^ 



Gn 



n -sr 



j: 2AD 

an* 

fi 



2ADn 



n ■sn-rr IM + 



(3n*-pn^ - Gn*-n( r + r^A 



2AD 



and so on, where 



■™ " (1) 



2oj e ~ AH 



D»D 



(3f 



(39) 



H 



[(1 - a)b -Gf3}^- + (Xb - G/3) 



(aP - Gb)^ + (hp - Gb) D 



C- BM^B 



C 



— (br + f3s) -p, 



A A 
a(l-a)-G 2 , D, = Xh-G 2 . 



(40) 
(41) 

(42) 
(43) 



From the definitions of the new parameters following Eq. ([H]), one may show that 
the sub-determinants D» and D± satisfy Dn > 0, D± > for all allowed values of 
a, X, h, and G, thereby justifying our use of formula in the derivation of Eq. (B5I). 



For the integral (P7|), the polynomial P[g, fc] is given by P = 2n*-q$/ri*-n. Since 
= a and /? = A from Eqs. ®), it follows from Eqs. (||), (||), (0), (g|), and © 
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that 



J[P] 



n -n 



n -p- 



{aX - Gb)' 



+ 



2ADw 



n 2 "e~ AH 



(44) 



H 



(b + X-G) (p 2 + m 2 ) - (br + As) -p. 



(45) 



3.4 Parameter Integration for the Matrix Method 



Once J[P] is known for a particular integral, we must complete the parameter 
integrations in Eq. fl28|) . For the two-loop integrals listed in the Appendix, we begin 
these integrations by changing variables from a\, ■ ■ ■ to the applicable subset of 
the new parameters A, j3,G,b,h,a and A, defined after Eq. (|31"|). For the sample 
integral (p7|), we find from Eqs. (EO) that 



roo roo /"1/2 rl—G rl—a ra roo 

da ± da 2 ... J = dG da dX db JA^dA 

JO JO JO JG JG JG JO 



(46) 



The integration ranges of the "finite" parameters A, /3, G, b, h, and a depend on which 
F factors are present in the original integral. For instance, if the factor containing 
q + k in Eq. ( ^7) did not also contain p, the b integration in Eq. ( f46"D would run from 
to a — G, rather than from G to a. In the event of a repeated F factor, there will 
be an additional finite parameter which can be integrated out immediately, since J 
will not depend on it. 

We see from Eqs. ( f4"5D and ( fI5D that D\\, D±, and H are independent of A, so that 
the A integration in Eq. (|46D is straightforward. Applying formula ([13|) to Eq. (0), 
we obtain 



J[P] A*dA 



71 



2^ 



n -n 



where 




D 2 D l r e H l+2e \D\\ D 



r(l + 2e)J 1 + r(2e)J 



Jo 



e = 2-u, 



*aH 



-26 



(47) 
(48) 



Note that the use of formula (|13l) requires H > 0. We can see that this requirement 
is satisfied in general by examining the way in which H was constructed: since the 
F factors and a parameters in Eq. (|28"D are non-negative, we find that the exponents 
in (|29|) and (^) are < for all values of q and k. Consequently, the sum of the 
exponents in Eq. fl34|) is < for all values of y. Taking y = 0, and noting that A is 
just a sum of a parameters, we see that H in Eq. (|42| ) is > for all allowed values 
of the old or new parameters. Furthermore, since r M and s M are independent of the 
quark mass m, one can deduce from Eq. ([45]) that H = only if b + X — G = 
(m ^ 0). We emphasize that the relationship b + A — G = applies specifically to 
the integral (f2~7|). 
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The presence of the divergent factor T(2e) in Eq. ([17]) reflects the UV- divergence 
of integral ( p7[ ) with respect to q and k combined. Additional divergences, known 
as subdivergences, may be found when we complete the remaining integrations from 
Eq. (|l6|), because Jo and J\ go to infinity at the boundaries of the integration region 
where a — ► and a — > 1 (corresponding to \k\ — > oo and |g| — ► oo, respectively). In 
the current example, each of these boundaries has three fewer dimensions than the 
full finite-parameter space, because from the limits of integration in (|46|) it follows 
that 

a > b > G > and l-a>A>G>0. (49) 

Thus, near a — 1, for instance, we could transform dGdX to (1 — a) 2 dxdy (with 
x and y being /imte parameters such that A = x(l — a) and G — y(l — a)). A 
subdivergence can therefore occur at a = 1 only if the integrand J or J\ diverges 
there at least as fast as (1 — a) ~ 3 . 

From Eqs. (|43|) , fljgp , and (|49D , we find that D\\,D± — > linearly as a — > 1, while 
the numerator of Ji goes to zero quadratically, and if remains positive unless 6 — > 
as well (assuming m ^ 0, as discussed above). Thus, we find that J is of order 
(1 — a) e_3 near a = 1, b ^ 0, while Ji is of order (1 — a) e ~ 2 there. Hence, by the 
criterion given above, the integral of Jo over the finite parameters diverges as a — > 1 
(and e — > 0), while the integral of J\ does not. Similar analyses at other boundaries 
of the integration region (including the case a = 1, b = 0) show that there are no 
other subdivergences in this example. Since we are only interested in finding the 
divergent parts of I2, and since J\ gives no subdivergences and is not multiplied by 
a divergent Gamma function in Eq. (f4~7D, we may drop J±. 



Due to the factor T(2e) in Eq. (0), the subdivergence in the integral of Jo will 
contribute a double pole to I2, while the finite part of the same integral will con- 
tribute to the single pole. The factor H~ 2e in Eq. ([48D complicates the integration 
of Jo, but some simplification is possible because this factor affects the pole parts 
of I2 only by way of the subdivergence at a = 1. To begin the simplification, we 
rewrite H~ 2e in (^) as 1 + (7J~ 2e — 1), and then expand the factors, multiplying 
the term (H~ 2e — 1), in powers of 1 — a, A, and G, using the definitions for Dm and 
Dx from Eq. (|KJ). In this fashion we obtain 

J »-z^ + TrS^[ 1 + 0(1 - a ' A ' G) l' (50) 

By counting powers of 1 — a, A, and G, we see that the integral of the term involving 
(H~ 2e — 1)0(1 — a, A, G) has no subdivergence; to find its finite part we therefore 
can set e = in this term (which happens to make it vanish). Accordingly, we may 
replace Jo by Jq(H), 



MH) = DfDV + (1 - aW-' ' (51) 
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without affecting the divergent parts of 12- 
The next step is to replace H in Eq. fl51~l) by 

E x = MmH = b(p 2 + m 2 ) - b 2 p 2 , 



where the result on the right follows from Eqs. ((H), (pfOD , ([43|), (^9|) , and from the 
equality (3 = A in Eqs. ([32]). To justify this replacement, we use the same kind 
of reasoning that led to the inequality (|25"1); in the present case, an appropriate 
inequality is given by 

< \H~ 2e -#r 2e | < 2t\H-H l \ (H' 2 ^ 1 + H{ 2e - 1 ), 

which holds for all allowed values of the parameters. Exploiting this new inequality, 
and noting that \H — H\\ goes to zero as a — > 1, and that both H and H\ go to zero 
(linearly) only when b — > 0, we can show that the integral of Jq{H) — Jq(Hi) , over 
the finite parameters, has no subdivergence and vanishes when e — > 0. In summary, 
we can replace H by Hi in Eq. (|5l|) without affecting the divergent parts of 12- 

Making this replacement, and integrating over b in accordance with Eq. (0), we 
obtain 



G 



D^i- e (1-a) 2 



1 /-G 
JO 



^r 2e - l)d6. (52) 



As a — > 1, the ranges of the second and third integrals in square brackets shrink 
to points, so that these terms contribute no sub divergences to the total integral 
over the remaining parameters. Hence, in the derivation of the finite parts of the 
integrals of these terms, we are allowed to set e = in their integrand (causing it to 
vanish). To facilitate the remaining b integration in Eq. (|5*2]), we shall express H^ 26 
in the form (cf. Eq. ([Tj 



H^ e = [b(p 2 + m 2 )-b 2 p 2 }- 2€ = (p2 + m 2-bp 2 r 2€ + b- 2e -l + 0(e 2 ). 

It now remains to integrate the factor multiplying the b integral in Eq. (Fj^) over 
A, a, and G, in accordance with Eq. ((46|). But this task is easy: since the b integral 
is of order e, only the divergent part of the integral over A, a, and G is needed. 



Finally, we must also integrate the first term on the right-hand side of Eq. ([52D 
over A, a, and G, again in accordance with Eq. (^). Because of the factor T(2e) in 
Eq. fl47[), both the divergent and finite parts of this triple integral will be needed. We 
begin the integration process by defining the new variables x = D±/a 2 , y = D\\/a 2 , 
and v = G/a, so that 

/ 1/2 dG /'"da r<iX a -^- = fiv Ay I" ^if^d*. 
Jo Jg Jg D 2 D± Jo Jv-v* y Jv-v* y 2 x 1 ~ e 



18 



From the definitions of y, v, and Dm, we find that a 2e = (y + v 2 + l)~ 2e . But 
this factor can be simplified because it affects the divergent and finite parts of the 
integral only at subdivergences. The subdivergence at a = 1 now manifests itself at 
v — y — x — 0. Near this point, we have 

(y + v 2 + l)- 2t = 1 - 2e(y + v 2 ) + e(l + 2e)(y + v 2 ) 2 . (53) 

Only the first term of this series gives a subdivergence, so we may set e = in the 
other terms (causing them to vanish, as usual). Away from the subdivergence, the 
integral is finite, so we can take e = in that region as well. Thus, in the current 
example, the factor a 2e may be replaced by 1. The integrations over x, y, and v 
pose no further problems. 

For a subdivergence at a = 0, we would use 

1 - 2e(v 2 + l)/y + e(l + 2e)(v 2 + 1) 2 /y 2 -...], 

valid for large y. If the subdivergence comes only from the first term in this series, 
then a 2e may be replaced by y~ 2e . 



(y + v 2 + l)- 2€ = y~ 



3.5 Reduction of Subdivergences 



The subdivergence in integral (]27|) is of the mildest possible nature; that is, J goes 
to infinity only just fast enough to make the finite-parameter integral diverge as 
a — > 1 and e — >• 0. For this reason, we were able to drop the 0(1 — a) term in 
Eq. (|50|), replace H by Hi, drop the last two integrals in Eq. (|52|), and drop all but 
the first term of series (0). For an integral with a more severe subdivergence, it 
may be helpful to begin with a partial calculation by the nested method, thereby 
reducing the degree of divergence with respect to one of the loop momenta before 
applying the full matrix method. 



A good example is the Euclidean-space integral 

P^(q,n,p) k^k v d 2uj qd 2uj k 



[(q - kf + m 2 } [k 2 + m 2 ] F 3 (q)F 4 (q) . . . 



(54) 



in which P^ v and i^i^, ... do not depend on k. Because of the quadratic subdi- 
vergence as | A; | — > oo, a direct application of the matrix method would lead to a 
finite-parameter integral proportional to / a e ~ 2 da near a = (after all parameters 
except a had been integrated out). To avoid having to deal with such a severe sub- 
divergence, we proceed in the spirit of Eq. ( |TTD by parametrizing the fc-dependent 
denominator factors only, to get 

h = I §^ f'dv H AY[k?kr\dA, (55) 

J -T3-T4 ... JO JO 
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where 



Y[f(k)} = Jd 2uj kf{k)exp(-A[k 2 -2vq-k + vq 2 + m : 



(56) 



Since Eq. ( p6|) has the same form as Eq. (|30|) (with z — k, M = A, B = Aug, and 
C = A{vq 2 + m 2 )), application of Eqs. fl35f) and ( |36f ) yields 



7T 

A, 

2\^2 i ^2 



Y[k^k l 



v 2 q»q v + i^A^ 1 



-AL 



v )q + m . We next use formula 
n UJ T(e)E^ 



and Eq. (p6|) to obtain 



/*oo /*oo 

^ / ^y[i]dA = / AY[EndA, (57) 

JO JO 



with L = (? 

AF[Pr]dA 

where = v 2 q^q u + L5^/(2e - 2). 

From the left and right ends of Eq. fl57|) , we see that -Z3 is unchanged if we 
replace the fc-dependent argument k^k v in Eq. (^3j) by the parameter-dependent 
polynomial E^ v (q,v), thereby reducing the degree of the sub divergence. We then 
cancel F factors, wherever possible, with factors in the terms of P^E^ ', parametrize 
the remaining F factors in accordance with formula (|9|), and finally complete the 
integration by using either the matrix method or the nested method, whichever is 
easier to apply. 



4 Renormalization 



Wherever a quark line appears in a physical process, any one of the diagrams of 
Figure 1, as well as higher- loop diagrams, could also appear. Hence, the effective 
quark propagator iSapScR is the sum of the bare propagator i5 a pS from Subsection 
2.1, plus contributions from self-energy processes with one loop, two loops, and so 
on: 



iS. 



cff 



iS + iS^iiS + iSiY^iiSY.i + S 2 )iS + . . . 



iS. 



eff, 



(58) 
(59) 



where £1 is the one-loop amputated Green function defined by Eq. (^), S 2 is the 
sum of the Green functions for the two-loop diagrams of Figures la to le, and so on. 
(£2 will also include amplitudes of one- loop processes with counterterm vertices, as 
explained below.) Since £,,■ is proportional to g 2 f j = 1, 2, 3, . . ., the right-hand side 



of 



is a power series in g . Similar series may be constructed for the effective 



gluon propagator and vertex factors: 
*5 a6 (GWU 



-r-rabap , -rrafrcrp 
+ii (2) _ 



iG pv + 



(60) 



V 1 eft) a /3 



1 (1)00 



1 r«f 

+ 1 (2)a/3 



+ . . . 



(61) 
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etc., where II(i), 11(2), T(i), and T( 2 ) are the amputated Green functions of the one- 
loop processes shown in Figures 2b to 2e. For the light-cone gauge, we have, in 
Minkowski space || [IT], [lSj , 



^i(p) = -ia s C F ( 2m — $ + 



1 (1)00 



finite, 



n -n 



Nc 
2 



■^n* M — r/i* 



+ finite, 



n - n 



(62) 
(63) 



■fin** 1 + T^n* 1 n*-kifcn^ 



& ~2 



+ 4 



n -n 



ngf (g) = za^^c 



11 



n • n 



n-kn*-n J 
g 2 n^ 

g M - - — 

n-q 



+ finite, (64) 



n* g 
n-q 



+ 



? - 



g 2 n" 
n-g 



71*"- 



n -qn h 
n-q 



+ finite, 



(65) 



<f(g) 



-iasS^-Nf^g^ - q»q v ) + finite; 



(66) 



here, p, g, and are external momenta as shown in Figure 2, Cp = "color factor" 
= | (<5afcCV = T®qT% k ), Nf = number of quark flavors = 6, Nc = number of colors 
= 3 (5 ad N c = f abc f bcd ), and 

5 s = /r(2-a;)(47r)- w 
in the "modified minimal subtraction" scheme, MS. 
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a) ^(p) 



b) r (1) 



c) r (2) (k) 




22 



Most of the terms in Eqs. (|5q) , (|6~0D , and (|6l|) diverge as — > 2. To make the 
effective propagators and vertex factors finite, we modify the Lagrangian density (0) 
in such a way that the bare propagators and vertex factors S, G, ig^, etc. will be 
replaced by renormalized propagators and vertex factors. For instance, if we let 



m 



m — 5m 



and 



ip + Nip 



(67) 



in Eq. ([I]), where N and 5m are both of order g 2 , then the quark propagator in 
Eq. (||) becomes 



S iea (p) = [(1 + 7°A rt 7°) tf-m + 5m) (1 + iV) 



(6* 



S(p) - S(p) [5m + 7 °iVt 7 (^ - m) + - m)iV] 5(p) + 0(# 4 ), (69) 



where the id term has been suppressed for clarity. At the same time, we see from 
Eq. (|58|) that in order to make S e s finite, we need to replace S by 



St, 



S-tS 



divergent part of Si S + 0(g 4 ) + 



(70) 



Comparing Eqs. fl69[ ) and (|70| ), substituting for S x from Eq. (|62|), and noting that 
tP = t/i* and 7°^* = ^7°, we find that we can eliminate the one-loop divergence 



from the effective quark propagator by taking |LT 
5m = 3ma s CF + 0(g 4 ) and iV = 5 S CV 



Similarly, to eliminate two-loop divergences from S e s, we will need to compute £2, 
and so on. 

The substitutions m — > m—5m and ^ — > tp+Nip produce the necessary counter- 
terms in the Lagrangian density (H). From Eq. (ffT|) we see that the counterterms 
involving are noncovariant, but we may hide this noncovariance by absorbing the 
factor (1 + N) into the normalization of the quark field if). Similar renormalizations 
of the gluon fields A a ^ lead to cancellation of the noncovariant divergences in the 
effective gluon propagator and three- and four-gluon vertex corrections. Bassetto, 
Dalbosco, and Soldati (BDS) have pointed out some time ago that noncovariant 
divergences may be eliminated in this way from the effective light-cone propagators 
and vertex factors, at all orders of perturbation theory |TTJ . We shall verify this 
claim explicitly for the case of the effective quark propagator at two-loop order. 



4.1 One- loop Counterterm Subtractions 



The use of renormalized propagators and vertex factors in place of S, G 
etc. affects S e s through the E factors in Eq. (|58|), as well as through S directly. 
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because the £ factors involve amplitudes which depend on the propagators and 
vertex factors through the Feynman rules. In accordance with standard procedure, 
these "indirect" contributions to S c s are represented by subtraction diagrams, such 
as those depicted in Figure 3. The diagrams of Figure 3 are specifically constructed 
from those of Figs, la to le by collapsing the divergent one-loop subgraphs (shown 
in Fig. 2) to the counterterm vertices denoted by x in Figure 3. In each case, 
the counterterm vertex factor is minus the divergent part of the amplitude of the 
corresponding subgraph. These amplitudes are given in Eqs. (|62]) to (|66|). The terms 
proportional to or n u in Eqs. (|64p and (|65D may be dropped (see Figs. 3c and 3d), 
since they are orthogonal to the light-cone propagators of the gluons. 



a) 



b) 





(1) 




(1) 



c) 



d) 



n 





(2) 




e) 




(2) 



(2) 



Figure 3: One- loop subtraction diagrams. The vertex factors for "x" are minus 
the divergent parts of the indicated one-loop amplitudes E 1; fm, etc. 




Figure 4: One-particle-reducible quark self-energy diagrams. The vertex factor for 
x is given by [— divergent part of Si.] 



The one-loop amplitudes S 1; Tm, rV 2 ), il( 2 ), and Il( 1 ) are of order g 2 , so that the 
amplitudes of Figures 3a to 3e are all of order g A . These amplitudes belong, therefore, 
to £ 2 , along with the amplitudes of Figures la to le, of course. In principle, the 
amplitudes of the one-particle- reducible diagrams (Figure 4) also contribute to X 2 , 
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but it is easy to see that the divergent parts of these amplitudes cancel completely 
with one another, due to the factorizability of all four integrals. 

All of the diagrams of Figures 3 and la to le contribute nonlocal divergences 
to E 2 ; however, the nonlocal contributions from Figs, la to le are exactly cancelled 
by the nonlocal parts from the corresponding diagrams of Figure 3, as expected on 
theoretical grounds |§ 

4.2 Two-loop Counterterms 

The divergent parts of the amplitudes for Figures la to le and 3a to 3e are shown 
in Table 1, excluding the nonlocal divergent terms which cancel as noted above. To 
find the local divergent part of the amplitude for a particular figure, one multiplies 
the numbers in the applicable row of the table by the corresponding factors at the 
top of the table, then adds the resulting terms together, and multiplies the sum by 
the color factor at the right-hand end of the row. (We recall that e = 2 — u, and 
a s = g 2 r(e)(47r) -w for modified minimal subtraction.) The divergent part of £2 
is the sum of the results for all ten figures, as shown in the "total" section at the 
bottom of the table. 



23, M. 
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Table 1: Local divergent parts of the amplitudes for Figures la to le and 3a to 3e. 
Each number is to be multiplied by the expressions at the top of its column and 
right-hand end of its row. C F , N c , Nf, and a s are defined after Eq. (67). Columns 
with 1 at the top give double-pole terms; columns with e, single poles. Vectors are 
in Minkowski space, and e = 2 — uj. 

n*-pjfi n-pr/i* 
m p 1 n*-n n*-n 



Figure 


1 e 


1 e 


1 e 


1 e 




rainbow 


la 
3a 


-8 16 
16 -24 


-3 -2 

6 2 

6 2 


-6 -40 
12 24 


6 8 

-12 -8 


^ioc 2 C F 
^ioc 2 C F 


over- 
lapping 


lb 
3b 


4 8 
-8 


4 2 

-8 -4 


14 60 

-28 -40 


-6 -12 
12 8 


\ia 2 s {Cl-\N c C F ) 
\ial{Cl-\N c C F ) 


spider 


lc 

3c 


4 -16 
-8 


-2 
-4 


-6 -20 
12 24 


-2 4 
4 8 


\ia 2 s N c C F 
\ia 2 s N c C F 


gluon- 
bubble 


Id 
3d 


44 124 
3 9 

3 u 


10 47 

3 9 

20 44 
3 3 


16 284 
3 9 

32 104 
3 3 


32 292 
3 9 

64 88 
3 3 


\ia 2 N c C F 
\ia 2 s N c C F 


quark- 
bubble 


le 
3e 


8 40 
3 9 

-16 
3 u 


4 14 
3 9 

8 8 
3 3 


8 64 
3 9 

16 16 
3 3 


8 64 
3 9 

16 16 
3 3 


\ia 2 s N f C F 
\ia 2 s N f C F 


total < 

(S 2 ) 




8 40 
3 9 

44 340 

3 9 

4 


4 10 
3 9 

22 49 
3 9 

-1 -5 
2 


8 16 
3 9 

44 172 
3 9 

-8 4 


8 16 
3 9 

44 172 
3 9 

-4 


\ia 2 s N f C F 
\ia 2 s N c C F 
^i<y 2 Cp 



A knowledge of the divergent part of £2 enables us to verify explicitly that 
the effective quark propagator may indeed be rendered finite at two-loop order by 
means of the substitutions m — ► m — Sm and ip — > tp + Nip , as implied by Bassetto, 
Dalbosco and Soldati. The 0(5 S ) parts of iV and Sm are already shown in Eqs. (ffl"l). 
Multiplying Eq. (|5Up by S^ 1 from the left and by (S^) -1 from the right, then solving 
for (5 c ff)~ 1 and replacing S by S rcn , we get 

(Seff)" 1 = (S^)' 1 - ffi! - i£ 2 - . . .. (72) 

To keep S e g finite at two-loop order, we must ensure that the divergent parts of 
E 2 in Eq. (|72|) are cancelled by terms from (S^en) -1 ; the latter depends on iV and 
Sm, as seen from Eq. (|68|). (The divergent parts of Si cancel already, thanks to the 
particular choice of the 0(5 S ) parts of iV and Sm.) Inverting and expanding the 
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right-hand side of (|68[) , we obtain 



(Srenr 1 = S' 1 + 5m + (5m -m)(N* + N + N*n) + N*tfN + N*]/> + rfN, (73) 

where N* = 7°iV"'"7 . Comparison of the terms of this expansion with the expres- 
sions at the tops of the columns in Table 1 suggests that the 0(5;?) parts of iV and 
5m should be similar in form to the 0(5 S ) parts shown in Eqs. ([H]). Therefore, let 
us try the ansatz 

5m = 3ma s C F + ma 2 s C F W +0(5;?), (74) 

and 

and then see if we can find expressions for W, X, and Y that will cause all 0(5;?) 
terms in Eq. ([72]) to cancel. 



Substituting from Eqs. ( [74]) and ([F^) into Eq. (|73|) , and applying the light-cone 
condition n 2 = n* 2 = 0, we find that the noncovariant 0(5;?) part of (S^en) -1 reads 



a 2 C 



F 



rr-n rr-n 



(76) 



S e s can be finite only if there is some value of X that makes expression (^) match 
the noncovariant, divergent part of «S 2 . Such a match is possible, but only because 
the coefficients under rf-piji in the "total" part of Table 1 are the negatives of the 
corresponding coefficients under n-ptf,*, except for the last two double-pole coeffi- 
cients in the last row. These last double-pole coefficients give rise to a term which 
cancels with the last term in brackets in (|76|) . This term in (|76]) derives purely from 
the leading term of N, via the N*p'N term in Eq. ([73]), so its coefficient cannot be 
adjusted (unlike X) to make it cancel with some arbitrary term from iY>2- Conse- 
quently, the anti-symmetry between the coefficients of n*-piji and n-piji* in £2 must 
be broken in exactly the way that we have found it to be broken, otherwise the BDS 
scheme would have failed. 

To complete our derivation of iV and 5m, having substituted from Eqs. ( |74] ) and 
(|75|) into Eq. (|73|) , we then substitute from Eq. (|73|) and Table 1 into Eq. ([72]), and 
set the total of the 0(5;?) terms on the right to zero. In this way we find the following 
expressions for W, X, and Y: 

N f ( 2 13 \ N c (11 221 \ C F n 11 



Y = Tl-3 + rJ + xlT-irJ+^U + r)' < 79) 
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Thus we have confirmed the claim of Bassetto, Dalbosco, and Soldati for the two-loop 
quark self-energy function. It is also significant to realize that the mass counterterm 
5m in the light-cone gauge, Eq. fl74|) , is exactly the same as in the general linear 
covariant gauge p2| , at least up to two-loop order. This result is a reassuring 
reflection of the gauge symmetry of QCD |25| . 



5 Conclusion 



In this paper, we have evaluated the complete quark propagator to two-loop or- 
der, and demonstrated for the first time its renormalizability in the noncovariant 
light-cone gauge. The chief technical tool in this three-year project was the pow- 
erful matrix integration technique which had been exploited in our first paper [Q 
to evaluate the divergent part of the overlapping self-energy function (see Fig. lb 
in the present article). Now, three years later, we have finally succeeded in com- 
puting, without approximation, the divergent parts of the four remaining two-loop 
diagrams in the quark propagator, namely: the quark-bubble diagram (Fig. le), the 
rainbow diagram (Fig. la), the spider diagram (Fig. lc), and the gluon-bubble dia- 
gram (Fig. Id), along with their counterterm graphs, depicted in Figs. 3e, 3a, 3c and 
3d, respectively. 

Other technical tools, apart from the matrix method and a reliable computer 
algebra program, included dimensional regularization, the n*-prescription for the 
spurious poles of (q-n)' 1 {n 2 = 0), as well as a detailed analysis of the boundary 
singularities in five- and even six-dimensional parameter space. 

The divergent part of the total two-loop correction S 2 to the quark propagator 
can be read off from Table 1, and has the form 



/II Ar 2 Ar \ _ / 1. n*-pi£- n-pii* 



+ ia 2 e{C 2 F ( -^ + 2 



3 



n*-pifi — n-piji* 



n*-n 



N f°F (, n 5 , A 
+ ( 10m- + 4 



N C C F ( 49 , 

+ -^—^ -85m + — p 1 - 43 
9 \ 4 r 



n*-pifi — n-pyji* 



n*-n 
n*-pifi — n-prji* 



n*-n 



JO) 



with a. 



g 2 Y{e)^f- 2 . 



The expression (|80|) has several noteworthy features: 
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(i) The anti-symmetry between the coefficients of n*piji and n-piji* (proportional to 
a 2 s C F ) is identical to the anti-symmetry predicted by the Bassetto-Dalbosco- 
Soldati renormalization scheme [Tcfl . 

(ii) The divergent part of S 2 is a local function of the external momentum p, even 
off mass-shell, since all nonlocal divergent terms cancel exactly. Notice that 
£2 contains both covariant and noncovariant components. 

(hi) The structure of S 2 implies the quark mass counterterm 8m, 

5m = 3ma s C F + ma 2 s C F W + 0(5*) , (81) 

which is gauge- independent, as expected (W is given by Eq. flTTD ). It is both 
interesting, and re-assuring from a calculational point of view, that the mass 
counterterm in the noncovariant light-cone gauge is exactly the same as in the 



general linear covariant gauge [22] - at least to two-loop order 



(iv) The factor N for the renormalization of the quark field ip is likewise local, 
albeit gauge-dependent: 

N = a s C F {t± - I) + a]C F ( X^± - Y )+ 0(5*) , (82) 



ji*-n 2y \ n*-n 

with X and Y defined in Eqs. (|78|) and (ff9|), respectively. 

The cancellation of nonlocal divergences, the consistency of our results with the 
Bassetto-Dalbosco-Soldati renormalization scheme, and the agreement between our 
mass renormalization and the two-loop covariant result all serve to demonstrate the 
reliability of the matrix integration method, as well as the validity, at two-loop order, 
of the n*-prescription for the unphysical poles of the light-cone gauge propagator. 

Since physically meaningful predictions can only be made with a theory whose 
divergences have been subtracted consistently, our result will help to make the light- 
cone gauge a more useful tool for practical two-loop calculations in QCD and other 
non-Abelian theories. Complete two-loop renormalization will also allow the com- 
putation of certain three-loop quantities, by means of renormalization group im- 



provements [2t . 



Of course, the cancellation of the divergences in the effective quark propagator 
is only part of the overall renormalization picture. In order to complete the two- 
loop renormalization of light-cone QCD, we must also find the counterterms which 
eliminate the two-loop divergences from the effective gluon propagator and vertex 
factors. We should be able to find these new counterterms using the same methods 
that have worked in this paper for the renormalization of the quark propagator. 
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Appendix 



The tables in this Appendix show one possible way in which the integrals I a to 
I e from Subsection 2.2 may be broken down into simpler integrals. Specifically, if 
we regard the columns of the tables as Cartesian vectors, then J a , for example, is 
the sum over all five tables of the dot product of the column headed by I a (if any), 
times the integral over q and k of the "integrand" column (plus some finite integrals 
which have been discarded). Blank spaces in the tables denote zeros. For example, 



The symbols in the denominators of the integrands are to be decoded as follows: 

n^n-q, q^q 2 , Q — > [(p - q) 2 + m 2 ], x — > \(p - q - k) 2 + m 2 ], 
o^n-k, k^k 2 , K — > [(p — k) 2 + m 2 ], y^(q-k) 2 . 

Note that all 4- vectors in this Appendix are in Euclidean space. 

For each row of each table, the divergent part of the integral over q and k of the 
expression in the "integrand" column is 7T 2uJ T(2e)(m 2 + p 2 )~ 2e times the expression 
in the final column, where e = 2 — u, and 






■n, 



p n *. n 



n 



* 






t = 



frf-p | (b-f)p 2 
n*n 2n-p 




P\\ 




P± =P~P\\, 



u = 



(m 2 + p 2 )d 
2n-p 




v 




(the "dilogarithm"). 
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4 


4 


h 


integrand, R 


(m 2 + p 2 ) 2e _ ^ /■ 2w 2 
7r 2w r(2e) P 7 y 


-4 






m 2 ^/nqQx k 


2w^ 


2 


-4 
4 


4 


i/tn-p/nqQx o 
ifin-p/ nqQ K o 


4(1 + 6 + cf)^p - 4(i + u)ifi 


-2e 
2 


2e 
-2 - 2e 
-2e 


2e 
-2 
-2e 


jji /n xKk 
^ / n Qx k 
^ /n Q Kk 


(| + 4-46 + 2c)^ 

(| + 4)^-4^ 

(1 + 4-46)^-4*^ 




-4 + 4e 


-4 + 4e 


^q-k/nqQ Kk 


(| - c)^ p + (w - t)^ 


4-4e 






^p-/c/nqQx k 


~~ D^p + (* + u ~~ 


4 


-4 




^p-g/nq xKk 
t^p-q/n QxKk 


(| + 3 - 2b)i P 
(| + 3- 26)^ p 


8-8e 


-4 + 4e 
-4 + 4e 


-8 + 8e 


1/tn-k/nq xKk 
^n-A;/n QxKk 


-f^ + j& - f# 




A A r 




(m 2 + p 2 )I/tn-k 
nqQxKk 




-4e 
4-4e 


4e 
4-4e 
-8 + 8e 
4-4e 


4e 

-8 + 8e 
4-4e 


^n-k/nq xKk 
gn-A;/nqQxK 
$n-/c/nqQx k 
^n-k/nqQ Kk 


-M + (2? + 1 - c )^; 
(-^-l + ^+^+(^ + l¥; 

+ (i + i-c)$ 


-16 + 8e 


-4 + 4e 
20 -8e 
-4 + 4e 
4-4e 


24 -8e 

-4 + 4e 
4-4e 


^n-p/nq xKk 
1/kn-p/n QxKk 
^n-p/nqQx k 
^n-p/nqQ Kk 


-il 

-il 

+ 

dp 1 


4e 
-8 


4-4e 

-4 
4 + 8e 
4-4e 
-4 + 4e 


4-4e 
4 

-4 + 4e 


^n-p/nq xKk 
$n-p/n QxKk 
^n-p/nqQxK 
^n-p/nqQx k 
^n-p/nqQ Kk 


(1 + 3-26)$ 
(I + 3 -26)$ 

2^ + (| + 3)^; 
2 ^ + (i + 3)#; 
2^+ (| + 4 -2b)i; 
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llKtigldllQ, -ft 


(m 2 +p2)2e r 
7r*T(2e) R 
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8(l-e)jrf+ 
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1 /qQxK 


±4-3-25 


4e#> + 
8(e - l)m 


2(2 - e)f+ 
2(1 -e)m 


8U - £)f + 
8(1 -e)m 


1 /qQx k 


± + 3-25 




2(e - 2)^+ 
2(e - l)m 


-2(1 + e)i 
—2em 


l /qQ Kk 


| + 4 - Ab 


-4e 


A-Ae 


A-Ae 


i /q xKk 






-8 + 8e 




^ / QxKk 
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integrand, R 
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